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The final analysis of the Extensive Air Shower (EAS) maximum (Xmax) depth distribution derived from the 
data of Tunka-25 atmospheric Cherenkov light array in the energy range 3 • 10^^ — 3 • 10^^ eV is presented. 
The perspectives of Xmax studies with the new Cherenkov light array Tunka-133 of 1 km^ area, extending the 
measurements up to 10^* eV, are discussed. 



1. INTRODUCTION 

The study of primary mass composition in the 
energy range 10^^ — 10^* eV is of crucial impor- 
tance for the understanding of the origin of cosmic 
rays and of their propagation in the Galaxy. The 
change from light to heavier composition with 
growing energy marks the energy limit of cos- 
mic ray acceleration in galactic sources (SN rem- 
nants), and of the galactic containment. An op- 
posite change from heavy to light composition at 
higher energy would testify the transition from 
galactic to extragalactic sources. Both changes 



*This work is supported by Russian Federation Ministry 
of Science and Education (Contract No 02.518.11.7073), 
the Russian Foundation for Basic Research (grants 07- 
02-00904, 05-02-04010, 06-02-16526) and the German Re- 
search Foundation DFG(436 RUS 113/827/0-1). Corre- 
spondence to v-prosin@yandex.ru 



are expected in the energy range of interest in 
the present investigation. 

To study the mean composition we use the rela- 
tion between the logarithm of mass In A and the 
depth Xmax of the extensive air shower (EAS) 
maximum: < X^ax >(x< \nA > - which is 
well-known from electromagnetic cascade theory. 
Xmax is derived for every event from the steep- 
ness of the atmospheric Cherenkov light lateral 
distribution function (LDP). We extend here the 
simplest method of performing such analysis (i.e. 
of relating the average < Xmax > and < \nA > 
values) to the more correct one based on the anal- 
ysis of the whole Xmax distribution. 

We summarize the results of the analysis of the 
Tunka-25 experiment, based on 25 stations cov- 
ering an area of 0.1 km^, and therefore tuned to 
operate between 3 • 10^^ and 3 • 10^'' eV. Both 
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quoted methods show the beginning of composi- 
tion change from hght to heavy at energies above 
10^^ eV. To extend the study to higher energies, 
higher statistics and consequently larger sensitive 
areas and solid angles are needed. Such an array, 
Tunka-133 (1 km^ area) is now under construc- 
tion close to its predecessor Tunka-25. Last win- 
ter the first part of the new array operated for 
about 270 h during clean moonless nights. The 
first results and perspectives of the new array are 
presented. 

2. TUNKA-25 FINAL RESULTS 

The detailed description of the Tunka-25 ex- 
periment and the procedure of primary energy 
measurement and ^mas analysis from the LDF 
steepness is given in 1 . The methods developed 
for the analysis provided the relative accuracy of 
energy to be ~15% 
than 30 g/cm^. 
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Figure 1. Depth of maximum, X^axi vs. primary 
energy Eq (7632 points). 
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Figure 2. Depth of maximum, Xmaxi distribu- 
tion for 6.4 < logio{Eo/GeV) < 6.5. Line - ex- 
periment, filled area - simulation for the com- 
plex mass composition containing '^70% of light 
(p+He) and ^-^30% of heavy (CNO+iron) nuclei 



The experimental plot of the depth of the 
shower maximum vs. primary energy Eq is shown 
in fig. 1. Data points are reported for energy 
£^0 > 2.5 • 10^^ eV at zenith angles 9 > 12° and 
for £^0 > 5 ■ 10^5 eV at 6* > 25°. In fact, the anal- 
ysis of possible distortions of the distribution has 
shown that no systematic errors are introduced 
in such energy-angular ranges. The distributions 
inside narrow logarithmic energy bins (0.1) have 
been analyzed. The analysis was done as follows. 
The experimental X„iax distribution is compared 
with the simulated one. The simulated distribu- 
tion is constructed from 4 partial distributions 
for different nuclei groups - p. He, CNO and Fe. 
Partial distributions are simulated with a " model 
of experiment" code assuming the QGSJET-01 
model of primary interaction. The code itself in- 
cludes all the essential parameter correlations and 
distributions extracted from CORSIKA and takes 
into account all the apparatus errors and selection 
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of events. A detailed description of this code is 
given in Ij . The weight of each group is selected 
for the best fit of the experimental distribution. 

The result for one of the logarithmic bins 
(6.4 < logw{Eo/GeV) < 6.5) is shown in fig. 2. 
We note that QGSJET-01 provides the best fit of 
the left edge of the distribution when compared 
with the other models. In principle the result of 
this procedure can give the relative weight of each 
group of nuclei within the total composition. But 
limited statistics and the relatively large width of 
the partial distribution make it almost impossible 
to distinguish 2 components inside the light (pro- 
ton+helium) and the heavy (CNO+iron) groups. 
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Figure 3. Mean logarithmic mass vs. primary en- 
ergy. Dense curve is the theory [2\ , dotted curves 
are the smoothing approximations of the experi- 
mental points. 



mental dependence of < In A > on primary en- 
ergy is shown in fig. 3, together with a theo- 
retical curve derived from [2]. 

One sees that the old method of interpolation 
shifts < InA > systematically by about 0.25 to- 
wards heavier composition, compared to the more 
strict method of the analysis of the full distribu- 
tion. The mean value obtained with the second 
method for the knee range of energies (3-10^^ eV) 
is close to that obtained in the recent balloon ex- 
periments for energy about 10 TeV [3]. The rise 
of < hi A > for energies above 10^^ eV is well 
visible. 

3. THE 1 km2 TUNKA-133 ARRAY 

To study the mass composition behavior in the 
energy range 10^^ — 10^^ eV, the new array Tunka- 
133 is under construction [4]. 
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A more stable solution can be obtained for 
the percentage of light (p-fHe) and heavy 
(CNO-|-iron) nuclei in the primary composition, 
that provides a stable estimation of the mean log- 
arithmic mass < In A >. The obtained experi- 



Figure 4. Plan of the Tunka-133 array. Black - 
4 clusters operated in winter 2007-2008, grey - 8 
clusters ready for operation in winter 2008-2009, 
open circles - the last 7 clusters to be deployed 
in 2009. 
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The array will consist of 133 detectors grouped 
into 19 clusters each composed of 7 detectors. 
The map of the array is shown in fig. 4. The 
new array provides much more information than 
the previous one. Each detector signal is digitized 
by an FADC with time step 5 ns. So the wave- 
form of every pulse is recorded, together with the 
preceding noise, as a total record of 5 /is duration. 

The minimal pulse FWHM is about 20 ns and 
the dynamic range of amplitude measurement 
about 10"*. The latter is achieved by means of 
two channels for each detector extracting the sig- 
nals from the anode and an intermediate dynode 
of the PMT with different additional amplifica- 
tion factors. 

Four clusters operated last winter between 
November and April. Data have been recorded 
over 270 hours during clean moonless nights. The 
average trigger rate was about 0.3 Hz, the number 
of the registered events was about 3 • 10^. 
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Figure 5. An example of an experimental event. 
The radii of the circles are proportional to the 
logarithm of the Cherenkov light flux 



4. RECONSTRUCTION OF EAS PA- 
RAMETERS. 

4.1. Structure of the data processing. 

The program of calibration and reconstruction 
of EAS parameters consists of three main blocks. 

1. The first block analyzes the primary data 
records for each Cherenkov light detector and de- 
rives three main parameters of the pulse: front de- 
lay at a level 0.25 of the maximum amplitude {ti), 
pulse area {Qi) and full width on half-maximum 
FWHM,. 

2. The second block of codes combines the 
data of different clusters and provides the rela- 
tive time and amplitude calibration. Data from 
various clusters are merged to one event, if the 
time difference for cluster triggers is less than 2 

flS. 

The time and amplitude calibration procedure 
is the same described in [1]. 

3. The third block of programs reconstructs 
the EAS core location, primary energy and depth 
of shower maximum. 



4.2. EAS core reconstruction Vifith the den- 
sity data Qi 

The first method of EAS core location recon- 
struction is based on the fitting of the Qi data by 
the lateral distribution function (LDF) with var- 
ied parameters of steepness (P) and light density 
at a core distance 175 m (Q175). This function 
was first suggested by the members of our collab- 
oration in [5 , and it is modified here to include 
large distance measurements: Q{R) — Qknf{R), 

rexp((«^(l-l-^)),i?<i?fe„ 
f{R)=l 200 >i?>i?..„(l) 

i (ff)"((2iT + l)/2)",i?>200 

Here R is the core distance (in meters), Rq is 
a parameter of the first branch of LDF, Rkn is 
the distance of the first change of LDF, Qkn is 
the light flux at the distance Rkn- The second 
change occurs at the core distance 200 m, b is the 
parameter of the third branch. This branch is 
checked till the distance 700 m with CORSIKA 
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simulated events. These 4 variables are strictly 
connected with two main parameters of the LDF 
- density at 175 m Q175 and steepness P: 

Qkn = Ql75(^fcn/175)^^-^ 

Ro = exp(6.79 - 0.564P), (m) 
i?fc„ = 207 - 24.5P, (to) 

4.84- 1.23Zn(6.5-P), P < 6 
3.43, P > 6 
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Figure 6. Lateral distribution resulting from fit- 
ting the measured light fluxes (points) with the 
expression (1) (curve) for the event from fig. 5. 



4.3. EAS core reconstruction with the 
pulse widths FWHM, 

In addition to the traditional method above de- 
scribed, a new method of EAS reconstruction us- 
ing Cherenkov light pulse FWHM has been de- 
signed and included into the analysis. To fit 
the experimental FWHMi, the empirical width- 
distance function (WDF) is used. It has a very 
simple analytic form (for FWHM > 20 ns and 
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Figure 7. Dependence of Cherenkov light pulse 
width on the core distance for the event from fig. 
5. Points - experiment, line - approximation by 
expression (2). 



R < 500 m): 

fl+lOO 

FWHM(P)^ll(™^7»°))^n. (2) 

FWHM(400) is related to the depth of EAS 
maximum Xmaxi that will therefore be recon- 
structed for each event by two independent meth- 
ods: from the LDF steepness P and the parame- 
ter FWHM(400). 

An example of a reconstructed shower is pre- 
sented in fig. 5. The LDF and WDF for this 
event are shown in fig. 6 and 7, respectively. 



PERSPECTIVES OF 
WIDTH ANALYSIS 



THE PULSE 



The absence of FWHM random fluctuations 
and a simpler expression for the WDF with re- 
spect to the LDF one seems to allow us applying 
the new method of EAS core reconstruction not 
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only inside, but also outside the array geometry, 
up to a certain distance. 

A similar idea of reconstruction of EAS core 
distance using pulse width was suggested many 
years ago by John Linsley [6]. But the use of 
the idea for charged particle detectors is prob- 
lematic because of essential random fluctuations 
of the signal form. Figure 7 shows that in case of 
Cherenkov light random fluctuations do not play 
an essential role. 

We realize that before using such method of 
core reconstruction, a detailed study of WDF up 
to core distances 1000 - 1500 m has to be per- 
formed not only with simulation but also experi- 
mentally. A positive result of this study would let 
us expand the sensitive area of the array by 5 - 10 
times compared with the geometrical area covered 
by the detectors for energies above 5 ■ lO^"^ eV. 

Such sensitive area increase will provide 20 - 
30 events with energies above 10"'^^ eV during one 
year (400 h during clear moonless nights) of ob- 
servation and ensure an overlapping of the Tunka- 
133 energy range with that of huge installations 
such as Auger. 
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6. CONCLUSIONS 

1. The analysis of Xj^ax distribution obtained 
with Tunka-25 array provides the mean logarith- 
mic mass of primary cosmic rays. For the energy 
range 3 ■ 10^*^ - 10^^ eV it is close to that ob- 
tained in balloon experiments for energies around 
10"'^^ eV and corresponds to a "light" composition 
^--^70% of p-l-He and ~30% of heavier nuclei. For 
energies above 10^® eV a rapid growth of < In A > 
is observed. 

2. The preliminary results obtained with the 
first stage of the new array Tunka-133 show 
unique possibilities of the new apparatus record- 
ing pulse waveform for each detector. This will 
provide more reliable measurement of X^ax with 
two methods in each individual event. 

3. The completion of Tunka-133 with its pro- 
jected area of 1 km^, together with the develop- 
ment of the new method of shower core recon- 
struction from pulse durations, will provide the 
reliable evaluation of < InA > up to 10^^ eV. 



